MOSER'S MATHEMAGICAL WORK ON THE EQUATION 

l k + 2 k + . . . + (m - l) k = m k 



PIETER MOREE 
In memory of Alf van der Poorten (1942-2010) 

Abstract. If the equation of the title has an integer solution with k > 2, then 
m > 10 10 . Leo Moser showed this in 1953 by amazingly elementary methods. 
With the hindsight of more than 50 years his proof can be somewhat simplified. 
We give a further proof showing that Moser's result can be derived from a von 
Staudt-Clausen type theorem. Based on more recent developments concerning this 
equation, we derive a new result using the divisibility properties of numbers in the 
sequence {2 2e+1 + 1}^L . In the final section we show that certain Erdos- Moser 
type equations arising in a recent paper of Kellner can be solved completely. 



1. Introduction 

In this paper we are interested in non-trivial solutions, that is solutions with k > 2, 
of the equation 

l k + 2 k + ... + {m-2) k + {m-l) k = m k . (1) 

The conjecture that such solutions do not exist was formulated around 1950 by Paul 
Erdos in a letter to Leo Moser. For k = 1, one has the solution 1+2 = 3 (and 
no further solutions). From now on we will assume that k > 2. Leo Moser [29] 
established the following theorem in 1953. 

Theorem 1. (Leo Moser, 1953). If (m,k) is a solution of ([T]), then m > 10 1C|6 . 

His result has since been improved. Butske et al. [6] have shown by comput- 
ing rather than estimating certain quantities in Moser's original proof that m > 
1.485 • io 9321155 . By proceeding along these lines this bound cannot be substan- 
tially improved. Butske et al. p. 411] expressed the hope that new insights will 
eventually make it possible to reach the benchmark 10 10 . 

The main purpose of this paper is to make Moser's remarkable proof of Theorem 
Q] better known, and with the hindsight and technological developments of more 
than 50 years, to give an even cleaner version of Moser's proof. This is contained in 
Section Moreover, we obtain the following refinement of Moser's result. 

Theorem 2. Suppose that (m, k) is a solution of ([I]) with k > 2, then 

1) m > 1.485 • 10 9321155 . 

2) k is even, m = 3 (mod 8), m = ±1 (mod 3); 
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3) m — 1, (m + l)/2, 2m — 1 and 2m + 1 are all square-free. 

4) If p divides at least one of the integers in (3), then p — l\k. 

5) The number (m 2 — l)(4m 2 — 1)/12 is square-free and has at least 4 990 906 prime 
factors. 

In fact, Moser proved (3) and (4) of Theorem [2] and weaker versions of parts (2) 
and (5). Readers interested in the shortest (currently known) proof of Theorem [2] 
are referred to Moree [25] . The deepest result used to prove Theorem [2] is Lemma 
[TJ Using a binomial identity due to Pascal (1654) a reproof of Lemma [1] was given 
recently by MacMillan and Sondow |18j . To wit, had Blaise Pascal's computing 
machine from 1642, the Pascalinejl worked like a modern computer, then Theorem 
[2] could have been already proved in 1654. 

In Section [3] we compare our alternative proof with Moser's original proof. 

In Section|4]we give a more systematic proof of Moser's result, which uses a variant 
of the von Staudt-Clausen theoremJl The relevance of this result for the study of 
the Erdos-Moser equation was first pointed out in 1996 by Moree (21] who used 
the result to show that the Moser approach can also be used to study the equation 
l fc + 2 fc + ... + (m-l) fc = am and a > 1 an integer. An improvement of the main 
result of [21] will be presented in Section [HJ 

The reader might wonder which other techniques have been brought to bear for the 
study of ([T]). Such techniques include Bernoulli numbers, considering the equation 
modulo prime powers, analysis (taking k to be a real, rather than an integer) and 
continued fraction methods. There is an extensive literature on the more general 
equation 

l k + . . . + (m - l) fc = y n , n>2, 

see, e.g., Bennett et al. [3]. That work incorporates several further techniques. 
However, those results do not appear to have any implications for the study of 
([T]). In Section [5J we give a taste of what can be done using Bernoulli numbers and 
considering ([T]) modulo prime powers. The main result here is Theorem 1 of [27] . We 
give a weakened (far less technical) version of this, namely Lemma HI Using that 
result and a heuristic assumption on the behavior of S r (a), a heuristic argument 
validating the Erdos-Moser conjecture can be given ([261 Section 6]). 

In Section [61 we consider implications for based on analytic methods, and 
in particular the recent work of Gallot, Moree and Zudilin [11] who obtained the 
benchmark 10 lt)7 and further improved this to 10 1C|9 by computing 3 ■ 10 9 digits of 
log 2. 

Section [7J is the most original part of the paper. Results on divisors of numbers 
of the form 2 2e+1 + 1 are used to show that if (m, k) is a solution of ([!)) such that 
m + 2 is only composed of primes p satisfying p = 5, 7 (mod 8), then m > 10 10 . 

In the final two sections we consider the Erdos-Moser variants 

l k + 2 k + . . . + (m - l) k = am k , respectively a(l k + 2 k + . . . + (m - l) k ) = m k 



The Pascaline was originally developed for tax collecting purposes! 
3 The proof given in Section 0] is implicit in Moree's [H] with a = 1. 
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(with a > 1 a fixed integer) and show that the latter equation (arising in a recent 
paper of Kellner [16]) can be solved completely for infinitely many integers a. 

This paper is in part scholarly and in part research. Leo Moser (1921-1970) was a 
mathematician of the problem solver type. For bibliographic information the reader 
is referred to the MacTutor History of Mathematics archive [30] or Wyman [40J . 

2. MOSER'S PROOF REVISITED 

Let S r (n) = Xlj=oi r - I n wna t follows we assume that 

S k (m) =m\ k > 2, (2) 

which corresponds to a non-trivial solution of ([I]). Throughout this note p will be 
used to indicate primes. 

Lemma 1. Let p be a prime. We have 

S r (p) = e r (p) (mod p), 

where 

I otherwise. 

Proof. Let g be a primitive root modulo p. In case p — 1 \ r we have 

^ . o'Cp-i) - 1 

Sr(p) = l^WY = : : — (modp), 

t-^ 1 9—1 

3=0 ^ 

and the numerator is divisible by p. In case p — l\r, we find by Fermat's Little 
Theorem that S r (p) = p — 1 = — 1 (mod p) as desired. □ 

Another proof using only Lagrange's theorem on roots of polynomials over Z/pZ 
can be given; see Moree [25] . The most elementary proof presently known is due to 
MacMillan and Sondow [18] and is based on Pascal's identity (1654), valid for n > 
and a > 2: 

£( n i 1 V fc(a)=a '* 1 - 1 - 

fc=0 ^ ' 

A further proof can be given using the polynomial identity 

p-i 

X V ~ Y - 1 = Y[(X - j) (mod p) 

3=1 

and Newton's identities expressing power sums in elementary symmetric polynomi- 
als. 

Lemma 2. In case p is an odd prime or in case p = 2 and r is even, we have 
S r (p x+1 ) = pS r (p x ) (mod p x+1 ). 
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Proof. Every < j < p x+1 can be uniquely written as j = ap x + /3 with < a < p 
and < (3 < p x . Hence we obtain by invoking the binomial theorem 

p— i p x — i p A -i p-i p x —i 

s r (p x+l ) = E E ^ + w = p E r + r ^ A E a E P_1 ( mod ^ 2A )- 

«=0 /3=0 /3=0 o=0 /3=0 

Since the first sum equals S r (p x ), and 2 X^a=o a = P(p ~ 1) = (mod p), the result 
follows. □ 

Proof of Theorem® Suppose that p\m — 1, then using Lemma [1] we infer that 

(m-l)/p-l p _ _ 

S k (m) = X>' + V )* = !?? —^(p) = ( m od p). (3) 

i=o j=i ^ ^ 

On the other hand m = 1 (mod p), so that by ([2]) we must have 

m — 1 



£fc(p) = 1 ( mod P)- ( 4 ) 
p 

Hence efc(p) (mod p), so that from the definition of e fc (p) it follows that efc(p) = 
— 1, and 

p\m — 1 implies p — l\k. (5) 
Thus (jlj) can be put in the form 

m — 1 

— — + 1 = (modp), (6) 

or 

m — 1 = —p (mod p 2 ). (7) 

We claim that m — 1 must have an odd prime divisor p, and that hence by (jSJ), 
must be even. It is easy to see that m — 1 > 2. If m — 1 does not have an odd 
prime divisor, then m — 1 = 2 e for some e > 2. However, by (J7J) we see that m — 1 is 
square-free. This contradiction shows that m — 1 has indeed an odd prime factor p. 

We now multiply together all congruences of the type ([6]), that is one for each 
prime p dividing m — 1. Since m — 1 is square- free, the resulting modulus is m — 1. 
Furthermore, products containing two or more distinct prime factors of the form 
(m — l)/p will be divisible by m — 1. Thus we obtain 



(m — 1) — h 1 = (mod m — 1), 



p\m— 1 



or 



V - + — - ee (mod 1). (9) 
^ v m - 1 



p m 

p\m— 1 



We proceed to develop three more congruences, similar to (Q, which when combined 
with OU) lead to the proof of part 1. Equation fl2]) can be written in the form 

S k (m + 2) = 2m k + (m + l) k . (10) 



THE ERDOS-MOSER EQUATION REVISITED 



5 



Using Lemma [T] and the fact that k is even, we obtain as before 

p\m + 1 implies p — l\k, (11) 



and 

m + 1 



P 



2 = (mod p). (12) 



From (I12p it follows that no odd prime appears with exponent greater than one in 
the prime factorization of m + 1. The prime 2 (according to H. Zassenhaus 'the 
oddest of primes'), requires special attention. If we inspect ([T|) with modulus 4 
and use the fact that k is even, then we find that m + 1 = 1 or 4 (mod 8). Now 
let us assume that we are in the first case, and we let 2^||m (that is 2^|m and 
2/+ 1 | m ). Note that / > 3. By an argument similar to that given in ([3]) we infer 
that S k (m + 1) = ff5'fc(2 / ) (mod 2 / ). Using Lemma El we see that S k (m + 1) = 
§S k (2 f ) = 2 /_1 (mod 2'), contradicting S fe (m + 1) = 2m k = (mod 2'). Thus 
m + 1 contains 2 exactly to the second power and hence (fT2|) can be put in the form 



777 -4- 1 

+ 1 = (mod p). (13) 

We multiply together all congruences of type ({TBI . The modulus then becomes 
(m + l)/2. Further, any term involving two or more distinct factors will be 

divisible by ^rj-^, so that on simplification we obtain 

V - + — — = (mod 1). (14) 
^ p m + 1 v ' K ' 

p\in+l 

We proceed to find two similar equations to ( !T4l) . Suppose that p\2m — 1, and let 
t = \( 2m ~ 1 — 1)- Clearly t is an integer, and m — 1 — tp + We have a k = (— a) k 
since k is even so that 2Sk( E Y-) = ^fc(p) (mod p) and hence, by Lemma [lj 

) = — r— (mod p). 



It follows that 



t-i p-1 (p-i)/2 



s'fcH = ^ + + z ' fc = (* + 2 )efc(p) (mod p) - (15) 

i=0 j=l i=l 

On the other hand 1 = (2m — 1 + l) k = (2m) k (mod p), hence m k ^ (mod p), 
so that (j2J) and ( fl5l) imply e^(p) 7^ 0. Hence p — l|/c, and by Fermat's little theorem 
m k = 1 (mod p). Thus (j2J) and ( fT5j) yield — (£ + |) = 1 (mod p). Replacing £ by its 
value and simplifying we obtain 

2m - 1 

+ 2 = 0(modp). (16) 

p 



6 



PIETER MOREE 



Since 2m — 1 is odd, (TTB|) implies that 2m — 1 is square- free. Multiplying congruences 
of the type (1161) . one for each of the r prime divisors of 2m — 1, yields 

2 r ~ 1 ((2m-l) - + 2) =0 (mod 2m- 1). 



V 

p\2m—l 

Since the modulus 2m — 1 is odd, this gives 

2 



p 2m 

p\2m-l 



(mod 1). (17) 



Finally we obtain a corresponding congruence for primes p dividing 2m + 1, namely 
(1191) below. For this purpose we write (T2]) in the form 

S k (m + 1) = 2m k . (18) 

Suppose p 1 2m +1. Set v = |( 2 "^ +1 — 1). Clearly v is an integer. We have m = pv+^- 
and find Skijn + 1) = {y + \)ek{p) (mod p). From this and ( Tl8|) it is easy to infer 
that ekip) = — 1, and so v + \ = — 2 (mod p). We conclude that 

p\2m + 1 implies p — l\k. 

Replacing v by its value and simplifying, we obtain 

h 4 = U (mod pj. 

p 

Note that this implies that 2m + 1 is square-free. Reasoning as before we obtain 

(mod 1). (19) 



p 2m + 1 

p\2m+l 



If we now add the left hand sides of (jUJ), (fT4|) . ( 1T7|) and ( TT9|) . we get an integer, at 
least 4. By an argument similar to that showing 2 { m, we show that 3 \ m (but in 
this case we use Lemma [2] with p = 3 and 3 A ||m and the fact that k must be even). 
No prime p > 3 can divide more than one of the integers m — 1, m + 1, 2m — 1 
and 2m + 1. Further, since m = 3 (mod 8) and 3 \ m, 2 and 3 divide precisely two 
of these integers. We infer that M\ = (m — l)(m + l)(2m — l)(2m + 1)/12 is a 
square-free integer. We deduce that 

1 2 2 4 111,, 

>4_ =3 (20) 



E 1 



p m — 1 m + 1 2m — 1 2m +1 2 3 6 

p\Nh 

One checks that (fT7|) has no solutions with m < 1 000. Thus (120]) yields (with 
a = 3.16) X] P |Mi p > a - F rom this it follows that if 



< a, (21) 

then m 4 /3 > Mi > Hp^^ anc l hence 



m > &*e B W\ (22) 
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with 6(x) = J2 P <x l°g£>, the Chebyshev ^-function. Since for example ( 12~TT) is satisfied 
with x = 1 000, we find that m > 10 103 and infer from (1201 that we can take 
a = 3| — 10~ 100 in fpZTT) . Next one computes (using a computer algebra package, say 
PARI) the largest prime pk such that J2 P < Pk \ < 3§, with pi,P2, ■ ■ ■ the consecutive 
primes. Here one finds that k = 4 990 906 and 

4 990 906 j . 

V — = 3.166 666 658 810 172 8584 < 3- - 10~ 9 . 

tr » 6 

This completes the proof of part 1 of the theorem; the remaining parts of the theo- 
rem have been proven along the way. □ 

Remark 1. Since for a solution of ((D), (m 2 - l)(4m 2 - 1)/12 has at least 4990 906 
distinct prime factors, it is perhaps reasonable to expect that each of the factors 
m — 1, m + 1, 2m — 1 and 2m + 1 must have many distinct prime factors. Brenton 
and Vasiliu [5] , using the bound given in part 1 of Theorem (2j showed that m — 1 
has at least 26 prime factors. Gallot et al. [H] increased this, using Theorem [5j to 33. 

Remark 2. Moser considered (EQ) modulo m — 1, m + 1, 2m — 1 and 2m + 1. Sondow 
and MacMillan |38] considered the equation also modulo (m — l) 2 and obtained some 
further information (this involves the Fermat quotient). 

3. Comparison of the proof with Moser's 

In this section we compare and contrast the proof of Theorem [2] with Moser's 
proof of Theorem [U 

Moser used only Lemma (U not Lemma [2J Consequently, he concluded that either 
m = 3 (mod 8) or m = (mod 8). In the first case we followed his proof but in the 
second case one has to note that we cannot use (j!4p . Letting M 2 = (m — l)(2m — 
l)(2m + 1) we get from ©, (PH), (US} 

- + + o 2 ; + 4 : > 3 - \ (23) 

p m — 1 2m — 1 2m + 1 3 

However, since 2 j M 2 , ( 1231) is actually a stronger condition on m than is ([2~0]) . 

The idea to use 3 f m, leading to a slight improvement for the bound on m, is 
taken from Butske et al. [Ej and not present in Moser's proof. (Actually they con- 
sider the cases 3 \ m and 3|m separately. We show that only 3 \ m can occur.) 

By using some prime number estimates from Rosser, Moser deduces that (12 ip 
holds with x = 10 7 and a = 3.16. In his argument he claims that by direct compu- 
tation one sees that (121]) holds with x = 1000 and a = 2.18. This is not true (as 
pointed out to me by Buciumas and Havarneanu). However, replacing 2.18 by 2.2 
in Moser's equation (21) one sees that his proof still remains valid. The present day 
possibilities of computers allow us to proceed by direct computation, rather than to 
resort to prime number estimates as Moser was forced to do. 

The advantage of the proof given in Section [2] is that it shows, in contrast to 
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Moser's proof and Butske et al.'s variation thereof, that every non-trivial solution 
satisfies the crucial inequality (120]) . 

4. A SECOND PROOF USING A VON STAUDT-CLAUSEN TYPE THEOREM 

In this section we show that Moser's four formulas (Q, (JHJ), (TTT1) and ffl~9]) can be 
easily derived from the following theorem. Indeed, using it a fifth formula can be 
derived, namely (126"]) below. 

Theorem 3. (Carlitz-von Staudt, 1961). Letr,y be positive integers. Then 

S (y) = V f = J° ( m ° d if r is odd; 

rV ~t 3 \-£ P -i|r, p\ v i ( mod V) otherwise. 

Carlitz [7] gave a proof of Theorem [3] using finite differences and stated that the 
result is due to von Staudt. In the case r is odd, he claims that S r (y)/y is an inte- 
ger, which is not always true (it is true though that 2S r (y)/y is always an integer). 
The author [20] gave a proof of a generalization to sums of powers in arithmetic 
progression using the theory of primitive roots. Kellner [15] gave a reproof for even 
r only) using Stirling numbers of the second kind. For the easiest proof known and 
some further applications of the Carlitz-von Staudt theorem, we refer the reader to 
Moree [25]. 



Second proof of Theorem® We will apply Theorem |3] with r = k. 

In case k is odd, we find by combining (I24p (with y = m) with ([I]) and using 
the coprimality of m and m — 1 that m = 2 or m = 3, but these cases are easily 
excluded. Therefore k must be even. 

Take y — m — 1. Then, using ([T]), the left hand side of (124"]) simplifies to 

S k (m - 1) = l k + 2 k + . . . + (m - 2) k = m k - (m - l) fc = 1 (mod m - 1). 
We get from ([24]) that 



E 



; +lE0(modm-l). (25) 



P 



p\m— 1, p — 1 1 A; 

Suppose there exists p|m — 1 such that p—l\k. Reducing both sides modulo p, we 
get 1 = (mod p). This contradiction shows that in (|2"5]) the condition p — l\k can 
be dropped, and thus we obtain ([S]). From ([S]) we see that m — 1 must be square- free 
and also we obtain ([9]). 

Take y = m. Then using ([1]) and 2\k we infer from (1241) that 

- = (mod 1). (26) 

p— l\k, p\m ^ 

Since a sum of reciprocals of distinct primes can never be a positive integer, we infer 
that the sum in (1261) equals zero and hence conclude that if p — l\k, then p\m. We 
conclude for example that (6,m) = 1. Now on considering (CD) with modulus 4 we 
see that m = 3 (mod 8). 
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Take y — m + 1. Then using (0Q) and the fact that k is even, the left hand side of 
241 simplifies to 



S k (m + 1) = S k (m) + m k = 2m k = 2 (mod m + 1). 

We obtain 

V ( m + 1 ) + 2 = (modm + 1), 

p\m+l, p— l\k 

and by reasoning as in the case y = m — 1, it is seen that p|m + 1 implies p — l\k, 
and thus (I14p is obtained. From Q14p and m = 3 (mod 8), we derive that (m + l)/2 
is square- free. 

Take y = 2m — 1. On noting that 

m— 1 

S fe (2m - 1) = J^(j fc + (2m - 1 - j') fc ) = 2S k (m) = 2m k (mod 2m - 1), 
we find that 

E(2m — 1) h . . 

^ ^ + 2m^ = (mod 2m -1). (27) 

p|2m— 1, p— l\k ^ 

Since m and 2m — 1 are coprime, we infer that if p\2m — 1, then p — l\k, m k 
1 (mod p) and furthermore that 2m — 1 is square-free. It follows from the Chinese 
remainder theorem that 2m fc = 2 (mod 2m — 1), and hence from ( 127|) we obtain 

(ED- 

Take y = 2m + 1. Noting that 

m 

S k (2m + 1) = ^(3 k + (2m + 1 - j) k ) = 2S k (m + 1) = Am k (mod 2m + 1) 



and proceeding as in the case y = 2m — 1, we obtain f[T9|) and the square-freeness 
of 2m + 1. To finish the proof we proceed as in Section [2] just below ffT9]) . □ 

With some of the magic behind the four Moser identities revealed, the reader might 
be well tempted to derive further identities. A typical example would start from 

4 * _ i* _ 2 k - 3 k = - ( mod m - 4 )- ( 28 ) 

P -i\k P 

p\m — 4: 

For simplicity let us assume that m ee 2(mod 3). We have (6, m — 4) = 1. For this 
to lead to a further equation, we need the left hand side to be a constant modulo 
m — 4. If we could infer that p\m — 4 implies p — l\k, then the left hand side would 
equal —2 (mod m — 4), and we would be in business. (For the reader familiar with 
the Carmichael function A, this can be more compactly formulated as \{m — 4)|/c.) 
Unfortunately a problem is caused by the fact that the left hand side could be 
divisible by p. Thus all we seem to obtain is that if m ee 2 (mod 3), and A(m — 4)\k 
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or 4 fc — l k — 2 k — 3 k and m — 4 are coprime, then 

1 2 

> = (mod 1). 

^— ' p m - 4 

p\m— 4 

In Section [TJ we will see that if we replace m— 4 by m+2 we can do a little better, the 
reason being that in this case, 2 k+1 + 1 appears on the left hand side, and numbers 
of this form have only a restricted set of possible prime factors. 

5. Bernoulli numbers and a cascade process 

Recall that the Bernoulli numbers are defined by the power series 

t _ Bkt k 
e* - 1 ~~ ^ k\ 

k=0 

They are rational numbers and can be written as B^ = Uk/Vk, with (t/jt, Vjt) = 1. 
One has B Q = 1, B\ = — 1/2 and B 2 j+i = for j > 1. By the von Staudt-Clausen 
theorem we can take for k > 2 even Vu = Yi P -i\kP m The Kummer congruences state 
that if k and r are even and k = r ^ (mod p — 1), then B^jk = -B r /r (mod p). A 
prime p will be called regular if it does not divide any of the numerators £4 with fc 
even and < p — 3. Otherwise it is said to be irregular. The first few irregular primes 
are 37, 59, 67, 101,.... 

The power sum S r (n) can be expressed using Bernoulli numbers. One has, see 
e.g. m (2-1)], 




Voronoi in 1889, see, e.g., [33| Theorem 2.8]), proved that if k is even and > 2, 
then V r S r (n) = U^n (mod n 2 ). From this result we infer that for a solution (k,m) 
of ([1]) we must have m\Uk and thus in particular v p (Uk) > v p (m), where we put 
u p (m) — f if p*||m. By a more elaborate analysis Moree et al. [27] improved this 
to VplBk/k) > 2v p [m). It shows (by the von Staudt-Clausen theorem) that if p\m, 
then p — 1 \ k (a conclusion we already reached using identity (1261) ). Invoking the 
Kummer congruences we then obtain the following result. 

Lemma 3. Let (k,m) be a solution of (pQ) with k > 2 and even. If p\m, then p is 
irregular. 

Let us call a pair (r, p) with p a regular prime and 2 < r < p — 3 even, helpful if 
for every a = 1, . . . ,p — lwe have S r (a) ^ a r (mod p). 

Lemma 4. If(r,p) is a helpful pair, and (k,m) a solution of (JTJ w'£/i fc even, £/ien 
we /iave k ^ r (mod p — 1). 

Proof. Assume that k = r (mod p — 1). By the previous lemma we must have 
p \ m. Now write m = m^p + b. Thus 1 < b < p — 1. We have, modulo p, 
Sk(m) = S r (m) = moS r (p) + S r (b) = SV(6). Thus if flTJ is satisfied, we must have 
^(ft) = 6 r (mod p). By the definition of a helpful pair this is impossible. □ 
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Since 2\k, and (2,5) is a helpful pair, we infer that 4|fc. Since (2,7) and (4,7) are 
helpful pairs, it follows that 6\k. From 4|fc and the fact that (4,17) and (12,17) 
are helpful pairs, it follows that 8\k. We thus infer that 24\k. It turns out that 
this process can be continued to deduce that more and more small prime factors 
must divide k\ for a detailed account with many tables see [26]. Given an irregular 
prime p and 2 < r < p — 3 even, one would heuristically expect that it is helpful 
with probability (1 — \jp) v ~ x which tends to 1/e, assuming that the values S r (a) are 
randomly distributed modulo p\ this is supported by current numerical data. 

Moree et al. [27] : using good pairs (of which the helpful pairs are a special case), 
showed that N x := lcm(l, 2, . . . , 200) divides k. Kellner [H] showed in 2002 that 
also all primes 200 < p < 1000 divide k. Actually Moree et al. (27J, p. 814] proved 
a slightly stronger result which combined with Kellner's shows that | k with 

N 2 = 2 8 • 3 5 • 5 4 • 7 3 • ll 2 • 13 2 • 17 2 ■ 19 2 ■ 23 • • ■ 997 > 5.7462 • 10 427 . 

An heuristic argument can be given suggesting that if L v := lcm(l, 2, . . . , v) divides 
k, with tremendously high likelihood we can infer that L w divides k, where w is the 
smallest prime power not dividing L v . It suffices that v > 11. To deduce that k is 
divisible by 24 is delicate, but once one has L v \k, there is an explosion of further 
helpful pairs one can use to establish divisibility of k by a larger integer. To add 
the first prime power w not dividing L v , one needs a number of helpful pairs that 
is roughly linear in v, whereas an exponential number (in v) is available. However, 
the required computation time increases sharply with w. This heuristic argument is 
the most convincing known to the author in support of the Erdos-Moser conjecture; 
details may be found in the extended version of [27J, [26] . 

Given a fixed integer a, one can try the same approach to study the equation 
Sk{ m ) — am k . Again one sees that once one manages to infer for example that 
120|/c, one can show that there must be larger and larger divisors. For many a, 
however, this 'cascade process' does not seem to 'take off' and it remains unknown 
whether all solutions with k > 2, if any, satisfy 120 1 k for example. 

If n — YliP? denotes the canonical prime factorization of n, then Q(n) = XI e « 
is the total number of prime divisors of n. Urbanowicz [39J proved a result which 
implies that given an arbitrary t, there exists an integer m t such that if (k, m) is a 
solution of ([T]) with k > 2 and m > m t , then Q(k) > t. 

6. The analytical approach and continued fractions 

Comparing Sk{m) with the appropriate integrals, it is easy to see that the ratio 
k/m must be bounded. A more refined approach gives 



On equating the left-hand side to m k and using (1 — l/m) m = exp(— 1 + 0(m 
one concludes that as m — > oo, we have 

log2 + 0( 



m \/m 
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By a rather more delicate analysis Gallot et al. [TTJ obtain that for m > 10 9 one has 
k / 3 C \ 

— = log2 1-- ^ , where < C m < 0.004. 

m \ 2m m l J 

As a corollary this gives that if (A;, m) is a solution of ([I]) with k > 2 and even, 
then 2k /(2m — 3) is a convergent Pj/<7j of log 2 with j even. This approach was 
first explored in 1976 by Best and te Riele [I] in their attempt to solve a related 
conjecture of Erdos; see also Guy (T2J D7]. The main result of [UJ reads as follows, 
where given N > 1, we define 

V(N) = {p : p — 1 | N} U {p : 3 is a primitive root modulo p}. 

Theorem 4. Lei N > 1 be an arbitrary integer. Let 

lo S 2 _r i- 1 
2jy — L a 0) a ii a 2i ■ ■ ■ \ — o-o H j 

ai H 

a 2 H 

5e i/ie (regular) continued fraction of (log2)/(2iV) ; with pi/qi = [ 
z-i/i convergent. 

Suppose that the integer pair (m,k) with k > 2 satisfies (JT]) | fc. Lei 

j = j(A^) &e i/je smallest even integer such that: 

a) a i+1 > 180iV-2; 

b) fe,6) = l; 

c) u p (qj) = v p {?> p 1 — 1) + v p (N) + 1 /or a// primes p e V(N) dividing q,j. 
Then m > q-j/2. 

Note that if for some integer iV we could prove that if all continued fraction digits 
<2j satisfy a, < lOOiV say and N\k, then (JTj) would be resolved! However, for a 
generic number £ G [0, 1] that is not a rational, one can show that the sequence of 
ai is not bounded above. The Gauss-Kuz'min statistics make this more precise and 
assert that the probability that a given term in the continued fraction expansion of 
a generic £ is at least b, equals log 2 (l + 1/6). Thus for a sufficiently large N, one 
expects that j(N) is quite large. This, in combination with the exponential growth 
of qj then ensures a large lower bound for m. (The numbers (log 2)/2N are expected 
to be generic.) 

The conditions b and c are of lesser importance. It seems that condition b is 
satisfied with probabilty 1/2. In practice, sometimes condition a is satisfied, but 
not b or c, and this leads to a larger lower bound for m. Condition c is derived 
using the Moser method, namely by analyzing the equation 



2(3 fc - l)(m - l) fc 
2m -3 




p-l\k 



that a solution (m, k) of (JT]) must satisfy. 

We leave it as a challenge to experts in the metric theory of continued fractions to 
determine the expected value of qj(N) on replacing (log2)/(2iV) above by a generic 
number £. Gallot et al. expect that conditions a and b lead to E (log qj(N)(£)) ~ c±N 
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and, taking into account also condition c, E(\ogqj(N)(£)) ~ c 2 N \og^ N for some 
positive constants ci,c 2 and (3. 

Crucial in applying the result is a very good algorithm to determine log 2 with 
many decimals of accuracy. Indeed, it is a well-known result of Lochs, that if one 
knows a generic number up to n decimal digits, the one can accurately compute 
approximately 0.97n continued fraction digits. For example, knowing 1000 decimal 
digits of 7i allows one to compute 968 continued fraction digits. 

Applying Theorem H with iV = 2 8 • 3 5 • 5 3 or iV = 2 8 • 3 5 • 5 4 , and using that N\N 2 
and N 2 \k, Gallot et al. obtained the current world record: 

Theorem 5. If an integer pair (m,k) with k>2 satisfies (pQ) ; then 

m> 2.7139 . 10*667658416 >10 io9 

Gallot et al. argue that, assuming one can compute log 2 with arbitrary precision, 
applying Theorem H] with N = N 2 should give rise to m > 10 10 . 

Interestingly, the results obtained by invoking Bernoulli numbers ('arithmetic') 
and analysis seem to be completely unrelated ('the arithmetic does not feel the 
analysis'). This strongly suggests that the Erdos-Moser conjecture ought to be true. 

7. A NEW RESULT 

This section focuses on new research; familiarity with the theory of divisors of 
second order sequences is helpful. The reader is referred to Ballot [2] or Moree [23] 
for more introductory accounts. 

Let S be an infinite sequence of positive integers. We say that a prime p divides 
the sequence if it divides at least one of its terms. Here we will be interested in 
the sequence S 2 := {2 2e+1 + 1}^L - It can be shown that p > 2 divides S 2 iff 
ord 2 (p) = 2 (mod 4), with ord 9 (p) (with p \ g) the smallest positive integer t such 
that g l = 1 (mod p). The set of these primes is known to have natural density 7/24 
[22]. Furthermore, if ord 2 (p) = 2 (mod 4) then 

p\2 2e+l + 1 iff 2e = ?lMP± _ 1 ( mod rd 2 (p)). (30) 

In some coding theoretical work the sequence 5*2 and its variants play an important 
role, as in [8j [13] and similarly in the study of the Stufe of cyclotomic fields 0, [22] 
and the study of Fermat varieties [311 [37] . 

If m + 2 is coprime with S 2 , then from (133]) and 2\k we can infer a fifth identity 
of Moser type, fl3"2"|) . This then leads to m > 10 10 for such m. We now consider 
the situation in greater detail. 

Theorem 6. Let N = (mod 24) be an arbitrary integer. Suppose that (m, k) is a 
solution of ([T]) with 

k>2, N\k and m < lO 10 ", 

then m + 2 has a prime divisor p > 3 such that: 

1) (ord 2 (p),iV) = 2,- 

2) k = ^^-l (mod ord 2 (p)). 
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In case m = 2 (mod 3), we can replace 10 10 by 10 10 . In case N = N2 we have 
p > 2 099. 

We first prove a corollary. 

Corollary 1. Suppose every prime divisor p > 3 ofm + 2 satisfies p = 5, 7 (mod 8). 

m>< 1nlo n (31) 
I 10 11 3 I m + 2. 

Proof. Using the supplementary law of quadratic reciprocity, (|) = (— 1)(p 2 ~ 1 )/ 8 ) one 
sees that if p = 5, 7 (mod 8), then ord2(p) ^ 2 (mod 4). Thus condition 1 is not 
satisfied, as for it to be satisfied we must have ord2(p) = 2 (mod 4). □ 



Put 

P(N) = {p> 3: (or d 2 (p),N) = 2}. 

We will study the set P{N) in greater detail with the ultimate goal of studying 
the iV-good integers, that is the odd integers n having no prime divisors in P(N). 
Note that in the proof of Corollary (TJ we established that integers composed only of 
primes p = 5, 7 (mod 8) are iV-good (with 24|iV). 

Corollary 2. Let N = (mod 24) be an arbitrary integer. If (m,k) satisfies (CQ), 
N\k and m + 2 is N-good, then m satisfies inequality ( 1311 . 

If p is to be in P(N), then p = 1 (mod 8) or p = 3 (mod 8). In the latter 
case we have ord2(p) = 2 (mod 4). In the former case it is not necessarily so that 
ord 2 (p) = 2 (mod 4), and numerically there is a strong preponderance of primes 
p = 3 (mod 8) in P{N). Indeed, we have the following result. 

Lemma 5. The relative density of primes p = 1 (mod 8) satisfying orc? 2 (p) = 
2 (mod 4) within the set of primes p = 1 (mod 8) is 1/6. 

Proof. We have seen that if ord2(p) = 2 (mod 4), then p = 1,3 (mod 8). If 
p = 3 (mod 8), then ord2(p) = 2 (mod 4). From this, the fact that <5(ord2(p) = 
2 (mod 4)) = 7/24 and the prime number theorem for primes in arithmetic progres- 
sion, we infer that the density of primes p = 1 (mod 8) such that ord2(p) = 2 (mod 4) 
equals ^ — | = ^. The sought for relative density is then i/| = ~. □ 

Thus if p = 3 (mod 8), then ord 2 (p) = 2 (mod 4) and if p = 1 (mod 8), then in 
1/6-th of the cases we have ord 2 (p) = 2 (mod 4). 

A further observation concerning the set P(N) is related to Sophie Germain 
primes. A prime q such that 2q + 1 is a prime, is called a Sophie Germain prime. 
Let denote the largest prime factor of M. 

Lemma 6. Let N = (mod 24) be an arbitrary integer. If q is a Sophie Germain 
prime, q = 1 (mod 4) and q and N are coprime, then p = 2q + 1 G P(N). 
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Proof. The assumptions imply that (J) = —1 and since p > 3 we infer that 
ord 2 (p) = 2q. Since (ord 2 (p), N) = (2q,N) = 2 we are done. □ 

There are 42 primes p in P(N 2 ) not exceeding 10 000. Of those 7 primes p are such 
that (p — l)/2 is not Sophie Germain, the smallest one being 7699. However, the 
Sophie Germain primes have natural density zero, whereas as we shall see P{N) has 
positive natural density. 

Given a rational number g such that g ^ { — 1,0, 1}, the natural density S g (d) of 
the set of primes p such that the order of g (mod p) is divisible by d is known to 
exist and can be computed; see e.g. Moree [23]. Using inclusion and exclusion one 
then finds that the set P{N) has natural density 



S(N) = J2(^(2d)-5 2 (4d)) l i(d), 



d\N 

where Nq is the product of the odd prime divisors dividing N and \x denotes the 
Mobius function. By Moree [23l Theorem 2] we then find that, for odd d, 

^2d)-5 2 (Ad) = l A H-J^- v 

p\d 

and hence 

V 



d\N p\d ^ p\N v ^ 

where we used that a multiplicative function / satisfies 

£m<o/w) = II(w(p))- 

d\N p\N 

Taking N = N 2 one finds that 

5 ^) = YA II (l--^) - 0.043 578 833- •■ 

2<p<1000 y 

By a result of Wiertelak, quoted as Theorem 1 in Moree (23], we have 

i = (i-5(iv))-^ + o^( r ^), 
P <*,7?p(v) h % x lo s x 

where the implicit constant may depend on N. From this result and [TU| Proposition 
4], we then infer that asymptotically the number of integers n < x that are A^-good, 
Nq{x), satisfies 

N G (x) ~ c N x\og~ siN) X, 

where 

1-6(N) / / v 

Xn{P> 



lim TT ( 1 - - 

t. — inn J- J- \ rr\ 



T(l - 8{N)) ^J; \ Pj \ P 
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with Xn{p) — if p = 2 or p is in P(N) and 1 otherwise. (As usual T denotes the 
Gamma-function.) Taking N = N 2 , a computer calculation suggests that c^ 2 ~ 0.54. 

Now if we have a sequence of random integers rij growing roughly as e 13 ^ for 
some constant (3 > 0, the integer rij is iV-good with probability log -5 ^ rij « 
c n{Pj)~ ^- The expected number of iV-good rij with j < x is then approximately 

The result that 2k/ (2m — 3) is a convergent Pj/qj of log 2 with j even and the result 
of Levy [17] that for a generic £ G [0, 1] that is not a rational 



log ft (0 



7T 



2 



lim = « 1.18, 

j-^oo j 12 log 2 

leads us to expect that the sequence m,j of potential solutions (kj, rrij) to (TjQ) coming 
from this result, is of exponential growth. Thus of the potential solutions (rrij,kj) 
with j < x, one expects about x 1-5 ^ 2 - 1 , that is roughly x - 96 , to be iVVgood. For 
those fl3T|) holds with m = mj. Thus if there would be say 10 10 potential solutions 
with m < 10 1Ql1 , then one expects roughly 3 • 10 9 to be iV 2 -good, and those can be 
excluded by Corollary [2j 



Remark. Given positive integers a, b, c, d, the density of primes p = c (mod d) such 
that p\{a e + & 6 }^L is known; see Moree and Sury [28]. Since 52 = {2 • 4 e + 1} 
that result cannot be applied to establish Lemma [5j 



oo 

e=0> 



Proof of Theorem [6j The idea of the proof is to show that if for every prime divisor 
p > 3 of m + 2 at least one the conditions 1 or 2 is not satisfied, then the identity 

= (mod 1) (32) 



E 



p m + ^ 

p\m+2 

holds. Using this we then show that m is bigger than the bound in the theorem; 
this is a contradiction. As usual we make heavy use of the fact that k must be even. 
We start with the equation 

2 k+1 + 1 = - Y (mod m + 2 ), (33) 

z — ' p 

p— l\k, p\m+2 

found on noting that Sk(m + 2) = 2m k + (m + l) k = 2 k+1 + 1 (mod m + 2) and on 
invoking Theorem El Suppose that p\m + 2. The idea is to reduce fl3"3"|) modulo p 
(except if p = 3, then we reduce modulo 9). 

If p = 3, then using Q\k we see that 2 k+1 + 1 = 3 (mod 9), and we infer that 
3 2 ||m + 2, that is we must have m = 7, 16 (mod 27). Next assume p > 3. 

First assume that ord2(p) ^ 2 (mod 4). Then p does not divide 5*2. Thus the right 
hand side of ( 133]) is non-zero modulo p, and this implies that p — l\k and p 2 \ m + 2 
and hence 2 fe+1 + 1=3 (mod p). 
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Next assume that orc^Qo) = 2 (mod 4), and condition 1 is not satisfied. Then 
ord2(p) and N have an odd prime factor in common, and by ( 1301) (with e = k/2) we 
get a contradiction to the assumption N\k. 

Finally, assuming that condition 1 is satisfied but not condition 2, the right hand 
side of (133]) is non-zero modulo p, and the same conclusion as before holds. By the 
Chinese remainder theorem we then infer that 2 k+1 + 1 = 3 (mod m + 2), and hence 
from fl33l we see that (1321 holds. 

Put M 3 = (m 2 — l)(4m 2 — l)(m + 2). By part 2 of Theorem [2] we infer that 
amongst the numbers m — 1, m + 1, m + 2, 2m — 1, 2m + 1, no prime p > 7 occurs 
more than once as divisor, the prime 2 occurs precisely twice, the prime 3 at most 3 
times and the prime 5 at most two times. Using this, we obtain on adding Moser's 
equations (|9|), (fl4|h (jTTj) and (fl9|) to (|32l): 



E 1 



1 2 2 4 3 109 

+ + - + - + - > — , (34) 



where 



p m — 1 m + 1 2m — 1 2m + 1 m + 2 30 

p\M 3 ^ 



109 12 1 

= 5 = 3.633 333 333 333 • • 

30 2 3 5 



Using the estimate 

^ - < log log x + 0.2615 H %— for x > 1, 

~ p log x 

p<x ° 

due to Rosser and Schoenfeld [35], (3.20)], we find that ^ p <«l/p < 3.633 32 with 
/3 = 4.33 ■ 10 12 . From another paper by the same authors [36] we have 

\6(x) - x\ < — ^ — , x > 678 407. 
1 v ' 1 401ogx' ~ 

Hence 

log(4m 5 ) > log(iV 3 ) > log ]Jp = 6(0) > .999/3, 

P</3 

from which we infer that m > 10 10 . 

In case m = 2 (mod 3) there are precisely two of the five terms m — l,m + 
1, 2m — 1, 2m + 1 and m + 2 divisible by 3, and in (|34|) we can replace 109/30 by 
109/30 + 1/3 = 119/30 = 3.966 666 • • ■ . In that case we can take = 4.425 • 10 17 
and this leads to m > 10 101 . 

The smallest two primes in P(N 2 ) are 2 027 and 2 099. For p = 2 027 we can 
actually show that condition 2 is not satisfied. To this end we must show that 
k ^ 1012 (mod 2 026). Computation shows that (1012,6 079), (3 038,6 079) and 
(5 064, 6 079) are helpful pairs. By LemmaHit then follows that k ^ 1 012 (mod 2 026). 
The smallest prime that possibly satisfies both condition 1 and 2 is hence 2 099. □ 
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Remark 1. We leave it as an exercise to the reader to show that (1311) can be refined 
to 



m > < 



' 10 1020 
io 1016 
io 1014 



10 



io 1 



if 3 
if 3 
if 3 
if 3 



m 
m 
m 
m 



5 
5 
5 
5 



m + 2 
m + 2 
2 
2. 



m 
m 



(35) 



In the same vein one can show that if (m, k) satisfies (pQ), k > 2 and m = ±1 (mod 15) 
or m = ±1 (mod 21), then m > 10 lo2 °. If, e.g., m = 1 (mod 15), then the sum in 
the left hand side of (JHJ) exceeds 1, so must be at least two. We infer that (12"01) holds 

with 3.1666... replaced by 4.1666 This then leads to m > 10 10 . The remaining 

cases are similar (they all lead to fl20|) with 3.1666... replaced by 4.1666....). 

Remark 2. Using the methods from Bach et al. p], it should be possible to compute 
the largest (3 such that ^ p< ^ 1/p < 109/30, respectively 119/30 exactly. They 
found, e.g., that the prime po = 180 124 123 005 660 046 7 is the largest one such that 



8. The generalized Erdos-Moser conjecture 
The Erdos-Moser conjecture has the following generalization. 
Conjecture 1. There are no integer solutions (m,k,a) of 

l k + 2 k + . . . + (m - l) fc = am k 
with k > 2, m > 2 and a > 1. 



(36) 



In this direction the author proved in 1996 [21] that ( 1361) has no integer solutions 
(a, m, k) with k > 1 and m < max(10 10 , a ■ 10 22 ). With the hindsight of more than 
10 years this can be improved. 

Theorem 7. The equation fl36|) has no integer solutions (a, m, k) with 

k>2, m< max (l0 91 ° 6 , a ■ 10 28 ) . 

Proof. (In this proof references to propositions and lemmas are exclusively to those 
in [21].) The Moser method yields that 2\k and gives the following four inequalities 

a + 1 



p-l\k 

p\m— 1 



+ — 

p m 



> 1, 



£ - + 

L — ' p m + 1 

p— i|fe r 



> 1. 



p\ m + 1 



E 

p— 1 1 fe 

p|2m — 1 



1 



2a 



2m — 



->i, v - 

1 - ' v 



2(a 



p— 1 1 k 
p|2m+l 



p 2m + 1 



> 1. 



(37) 



(3* 



Since p|m implies p — 1 \ k (Proposition 9), we infer that (6, m) = 1. Using this we 
see that Mi = (m 2 — l)(4m 2 — 1)/12 is an even integer. Since no prime > 3 can 
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divide more than one of the numbers m — 1, m + 1 , 2m — 1 and 2m + 1 , and since 2 
and 3 divide two of these numbers, we find on adding the inequalities that 

El a a + 1 2a 2(a + 1) 11 1 
- + + + + '- > 4 = 3-. 
p m — 1 m + 1 2m — 1 2m +1 2 3 6 

p— l\k, p\M\ 

Using that a(k + 1) < m < (a + l)(k + 1) (Proposition 2), we see that in the latter 
equation the four terms involving a are bounded above by 6/ (k + 1). Since k > 10 22 
(Lemma 2), we can proceed as in the proof of Theorem [2] and find the same bound 
for m, namely m > 1.485 • 10 9321155 . 

Earlier it was shown that if k > 1, then k > 10 22 . To this end Proposition 6 with 
C = 3.16, s = 664 579 = 7r(10 7 ) and n the 200-th highly composite number was 
applied. Instead we apply it with C = 3| - 10~ 10 , s = 4 990 906 and n the 259-th 
composite number c 2 so (this has the property that the number of divisors of C259 < s, 
whereas the number of divisors of c 2 6o exceeds s). Since n = C259 > 5.583 4 • 10 27 it 
follows that k > 2n > 10 28 . Since m > a(k + 1), the proof is completed. □ 

Remark 1. The above proof shows that if ( 1361) has a solution with k > 2, m > 2 
and a > 1, then m must be odd. An easy reproof of this was given by MacMillan 
and Sondow 



Remark 2. The reader might wonder whether the method of Gallot et al. can be 
applied here as well to break the 10 l ° 7 barrier. For a fixed integer a this is possible 
if one manages to establish that N\k with iV large enough. Gallot et al. showed 
that 2k/ (2m — 2a — 1) is a convergent with even index of log(l + 1/a) for m large 
enough. For a given a this can be made effective. Establishing that N\k along the 
lines of Section [5] is not always possible (see the last paragraph of that section). 

Challenge: Reach the benchmark 10 1()7 in Theorem [71 

9. The Kellner-Erdos-Moser conjecture 

Kellner (16] conjectured that if k, m are positive integers with m > 3, the ratio 
Sk(m + 1)/Sk(m) is an integer iff (k,m) G {(1,3), (3,3)}. Noting that Sk(m + 1) = 
Sk(m) +m k one easily observes that this conjecture is equivalent with the following 
one. 

Conjecture 2. We have aSk(m) = m k iff (a, k,m) G {(1, 1, 3), (3, 3, 3)}. 

If this conjecture holds true, then obviously so does the Erdos-Moser conjecture. 

It is easy to deal with the case m = 3. Then we must have a(l + 2 k ) = 3 k , and 
hence a = 3 e for some e < k. It follows that 1 + 2 k = 3 k ' e . This Diophantine 
equation was already solved by the famous medieval astronomer Levi ben Gerson 
(1288-1344), alias Leo Hebraeus, who showed that 8 and 9 are the only consecutive 
integers in the sequence of powers of 2 and 3, see Ribenboim [3U pp. 124-125]. This 
leads to the solutions (e, k) G {(0, 1), (3, 1)} and hence (a, k,m) G {(1, 1, 3), (3, 3, 3)}. 
Next assume that m > 4 and k is odd. Then by Theorem[3]we find that mim — l)/2 
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divides m k , which is impossible. We infer that to establish Conjecture |2J it is enough 
to establish Conjecture [31 where 

A = {a > 1 : aSk(m) = m k has a solution with 2|fc, k > 2, m > 4}. 

Conjecture 3. The set A is empty. 

The next result shows that if a = 2 (mod 4) or a = 3, 6 (mod 9), then a ^ A. 

Theorem 8. Let k > 2 be even. Suppose that q\a is a prime such that q 2 \ a and 
q — l\k. Then aSk{m) ^ m k and hence a ^ A. 

Proof. Suppose that aSkijn) = m k . Let g e ||m. Note that e > 1. Using Theorem 
[3] we find that Skijn) = ^Sk(q e ) = — — (mod q e ). Now we consider the identity 

aSk(m) = m h modulo q e+1 and find — = m k = (mod q e+1 ), contradicting 

q e+1 \\am. It follows that aS^m) ^ m h . □ 

Note that if a ^ A, then the equation aSk(m) = m k can be solved completely. 
The author is not aware of earlier 'naturally' occurring Erdos-Moser type equations 
that can be solved completely. He expects that further values of a can be excluded 
and might come back to this in a future publication. 
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they are back full force in the title of [25]. Also whilst at Macquarie, thanks to 
questions by then visitor Patrick Sole, I got into the study of divisors of a k + b k , not 
realizing there is a connection with the Erdos-Moser equation (as the present article 
shows) . 

Had Alf learned that the present record for solutions for EM is based on a con- 
tinued fraction expansion (of log 2), I am sure he would have been pleased. 

I had a wonderful year in Australia and will be always grateful to Alf for having 
made that possible. 
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